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Abstract 

In this study, the hybrid Cramer-Rao bound (CRB) is developed for target localization, to establish the sensitivity 
of the estimation mean-square error (MSE) to the level of phase synchronization mismatch in coherent Multiple-Input 
Multiple-Output (MIMO) radar systems with widely distributed antennas. The lower bound on the MSE is derived 
for the joint estimation of the vector of unknown parameters, consisting of the target location and the mismatch 
of the allegedly known system parameters, i.e., phase offsets at the radars. Synchronization errors are modeled as 
being random and Gaussian. A closed-form expression for the hybrid CRB is derived for the case of orthogonal 
waveforms. The bound on the target localization MSE is expressed as the sum of two terms; the first represents the 
CRB with no phase mismatch, and the second captures the mismatch effect. The latter is shown to depend on the 
phase error variance, the number of mismatched transmitting and receiving sensors and the system geometry. For a 
given phase synchronization error variance, this expression offers the means to analyze the achievable localization 
accuracy. Alternatively, for a predetermined localization MSE target value, the derived expression may be used to 
determine the necessary phase synchronization level in the distributed system. 

Index Terms 

MIMO radars, Hybrid CRB, mismatch parameters, localization. 



I. Introduction 

Improvement in target parameter estimation capabilities is a primary advantage of MIMO radar systems [l]-[5], 
In particular, target localization with coherent MIMO radar systems, utilizing widely distributed antennas, offers 
significant advantages [5]. Typically, performance analysis of system parameter estimation problems is based on the 
derivation of the Cramer-Rao bound (CRB), which sets a lower bound on the estimation MSE for unbiased estimators 
[6]. Such an evaluation is provided in [2], [5] for coherent MIMO radar systems, demonstrating a localization 
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accuracy advantage, inversely proportional to the signal carrier frequency. In addition, a spatial advantage of the 
order of the product of the number of transmit and receive radars is also incorporated in the CRB. 

This performance gain comes with the challenge of attaining phase synchronization in a distributed system. 
Errors introduced to the system parameters by phase synchronization mismatch, will result in parameter estimation 
mean-square error (MSE) degradation and bias. In this work, the hybrid CRB (HCRB) is used to test the sensitivity 
of the target localization MSE to phase errors. The HCRB takes into account deterministic unknown parameters, 
such as the target location, as well as random parameters, phase calibration errors, in this case. This method has 
been applied to passive source localization [8], [9] for the problem of source bearing and range estimation with 
uncertainty in the sensors' locations or phase synchronization errors. 

In this work, the HCRB is derived for coherent MIMO radars, with phase synchronization errors. A closed-form 
expression for the HCRB for the target's location (x, y) is derived, providing the means to assess the effects of 
phase errors on the localization accuracy. The effect of the number of radars, their geometric layout, and the phase 
mismatch MSE is incorporated in the HCRB terms. 

The paper is organized as follows: brief theoretical background is provided in Section II. The system model is 
introduced in Section III, and the HCRB on the targets localization estimation errors is derived. Numerical examples 
are presented in Section IV. Finally, Section V concludes the paper. 

II. Background 

The hybrid CRB provides a low bound on the MSE of any unbiased estimator for an unknown parameter(s), 
where the parameters are partially deterministic and partially random [7]. Given a vector parameter 6 = [9 nr , 9 r ] T , 
where 9 nr stands for the nonrandom parameter vector and 9 r for a random parameter vector, its unbiased estimate 
9 satisfies the following inequality [7]: 

Ee nr ,e r {(£ - 0*) ft - ^) T J > [Jjj 1 (0 nr ,9 r )] .. , (1) 

where Jh (9) is the hybrid Fisher Information matrix (HFIM) expressed as 

J if {9nr, r ) = J D + Jp (2) 

The elements of the matrices J d and J p given by 

d\n P {v\e nr ,e r ) 



[3D] id -E 6rl g nr E T]9nrfir I \\, (3) 



and 



r T , p f dHnp(0r\0 nr ) \ 



where p (r\9 nr , r ) is the conditional, joint probability density function (pdf) of the observations and p {9 r \9 nr ) the 
conditional joint pdf of 9 r - The matrix J/) represents the contribution of the data and the matrix Jp represents the 
contribution of prior information. 
The HCRB matrix is defined as 

HCRB = [J h (9 nr , 9 r )]~ 1 . (4) 
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In cases in which the observation statistic is expressed in terms of p (r|K nr , K r ), and the relationship between the 
unknown parameters 6 nr , r and K nr , n T is given by Kj = fj(0), the chain rule, can be used to express Jh {@nr, ®r) 
in an alternative form [11]: 




where the elements of the matrix P are given by [P] 4 • = 



III. HCRB with Phase Mismatch 



In this section, the HCRB is developed for target localization. A point target is assumed with complex reflectivity 
$ = $Re + j$im, located in a two dimensional plane at coordinates X Q = (x ,y ). Consider a set of M 
transmitting stations and N receiving stations, widely distributed over a given geographical area, and time and phase 
synchronized. A set of orthogonal waveforms is transmitted, with the lowpass equivalents Sfc (i) , k = 1, . . . , M, 
and effective bandwidths (3 [10]. The signals are narrowband in the sense that for a carrier frequency of f c , the 
narrowband signal assumption implies (3 2 / f% <C 1 

In [5], perfect phase synchronization was assumed. In practice, synchronization errors exists, modeled here as zero 
mean Gaussian random variables with standard deviation a\ and denoted by A</> = [A0 tl , A</> t2 , A.<f>t M , A<f> ri , A<j) r2 , A(f> rN } T ', 
where A(f> tk and A(f> rt are phase errors at transmitting radar k and receiving radar i, respectively. The phase errors 
introduced by the different stations are assumed to be statistically independent. The vector of unknown parameters 
is defined by 



— [0 nr , r ] T , 



(6) 



where 9 nr = [x , y , ^Ro, #im] denotes the deterministic unknowns and 9 r = A0 T denotes the random unknowns. 

The estimation process is based on the signals observed at the receiving sensors. The signal received at sensor 
I is a superposition of the transmitted signals, reflected from the target, and given by: 



A I 




(7) 



fe=i 



where accounts for the phase information and has the value of r/^ = exp (— j27r/ c r^fe) exp (— j {A<p tk + A<f> rt )). 
The noise ni (t) is assumed to be circularly symmetric, zero-mean, complex Gaussian, spatially and temporally 
white with autocorrelation function a^S (t). The propagation time, r^k, is a sum of the time delays from station k 
to the target and from the target to station £, and may be expressed as 




(8) 



where c denotes the speed of light, (x tk ,y tk ) denotes the location of transmitting radar k and (x r e,y r i) denotes 
the location of receiving radar I. The following vector notation is introduced: r = [th,ti2, ■■■, T ik, tnm] T ■ 



The received signals are separated at the receiver by exploiting the orthogonality between the transmitted 
waveforms. The signal in (7) is defined as a function of the time of arrival, t^, the reflectivity value and 
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the phase mismatch A(f>. The vector of unknown parameters for the observations ri (t) is expressed as a function 
of the time delays r rather than a function of the unknown location (x ,y ) (as seen in (8)); i.e., the vector of 



. The following notation 



unknown parameters is denoted by k — [K nr , k] T , with n nr — [r T , $r , #i m ] and k t = 
is defined for later use: r = [n (t) , . . . , r N (t)f ', Q = MN, L = M + N. 

In order to derive the HFIM given in (2) and (3), the conditional joint pdf p(r\n) is required. For the signal model 
given in (7), the conditional joint pdf of the observations (time samples at multiple receive antennas) parametrized 
by the unknown parameters vector k, is then 



p (r\n) oc exp 



(9) 




M 



r e(t) - ^2$s k (t - r lk ) r]e k 



fe=i 



dt 



The observation is given as a function of k. Therefore, the matrix P, defined following (5), needs to be derived. 
The relation given in (8) is used, resulting in 





P = 



1J 2xQ 





L(2+L)x(2+L) 



(10) 



with 



D = — - 

c 



-i T 



cos ai + cos 72 sin 012 + sin 72 



(11) 



cos aM + cos 7at sin ccm + sin -f^ 
where is the bearing angle of the transmitting sensor k to the target, measured with respect to the x axis, and 
-fi is the bearing angle of the receiving radar I to the target, measured with respect to the x axis. 

Using the conditional pdf p(r\n) in (9) and the Gaussian distribution of the phase errors, the HFIM J# (k), 
defined by (2) and (3), is derived in Appendix I, resulting in 



J/7 (^nn K-r) 



G H 



where matrices G and H are defined by 



G = 



JQx(2+L) 



and 



H 



(12) 



(13) 



(14) 



(2+L)x(2+L) 

and the other submatrices in (12), (13) and (??) are defined and derived in Appendix I (see (27), (29), (30) and 
(31)). Applying (10) and (12) in (5) yields 



3 H {9) 



DR T D T DG 
G T D T H 



(15) 
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The HCRB for the unknown parameters (x ,y ) may be derived from (15), applying the relation given in (4) [5]: 

HCRB (x , Vo) = [DR r D T - DGH X G T D T ] . (16) 
To find the closed-form solution to HCRB (x Q , y a ), the matrix H 1 is expressed using the formula for the inverse 



of a partitioned matrix [12]: 



H" 



i it 



[£,9 - F^aA a 1 FJ a ] 



[ H_1 ] 2 2 = -[ F ^A^ 1 Fj A -A A ]- 1 



(17) 
(18) 



and 



[H- 1 ] 12 = [H- 1 ] T 21 = E^F^ [F^S-^^a -A A ]^ (19) 

where Aa = ^Sa + ■ The term [E# — F^ A A A 1 F^ A ] 1 in (17), is transformed based on the formula for the 
inverse of a matrix B of the form B = A + XRY, given in [12]. Following some additional matrix manipulations, 
the HCRB for the location MSE can be expressed as 

HCRB (x ,y ) = [j F - JfPa'Jf] (20) 

= CRB (i„,|/ ) + ACRB, 

where CRB Q (x a , y a ) = J^ 1 is the CRB with no phase mismatch, and ACRB = \jp — J fPa 1 J f] 1 represents 
the increment in the bound due to phase synchronization errors. The matrices Jf and Pa are defined by 



J F = DR r D J DFrfE^F^D 



and 



Pa — DFt^S^ A R A 1 Fj A S l} 

- 2Rc{DF TA R A l Fj A S- l F T 1? D T } 



(21) 



+ df tA r a 1 f t a d , 

and the matrix R A X can be calculated recursively using the formula for the inverse of the sum of matrices [13], 



resulting in 



R A X 



N\1 n T 



AiImxM + M(l-JVAi) 





\ T i MA 2 1 1 T 

Mi-NxN + N{1 -MX 2 ) 11 



where 1 = [1, 1, , 1] and the terms Ai and A2 are 



Ai = 

y 2 snr crj 



and 



A 2 = ( 



M + „ 1 3 

2 snr ai 



Calculating the explicit value of ACRB, we get 



ACRB 



Jf — Jf I Mo ^ km'Bm J Jf 



m— 1 



(22) 



(23) 



(24) 
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where the constants fj, m , m = 0, 3 are functions of the phase synchronization error variance a\ (through Ai and 
A2, defined in (23)) and the number of transmitting and receiving radars M and N, as follows: 

8n 2 (f 2 + (3 2 ) snr 



Mo = 



(25) 



/xi = Ai/M + \ 2 /N, 
H 2 = Ai (N 2 /M) and fc 3 = A 2 (M 2 /N) . 
The matrices B TO , m — 1,2,3 depend on the geographical layout of the radars with respect to the target location: 



B, 



using the following D m matrices: 



\D m l]t,i 
[D ro l] M [D ro l] 21 



Di = cD 



[D ro l] M [D m l] 2jl 



(26) 



D 2 



and 



cos ai sma2 



cos olm sin cum 



cos 72 sin 72 



-, T 



J Mx2 
T 



Nx2 



D 3 

cos 7 iv sin 7tv 

The expression for the HCRB as given in (20), offers an interesting observation on the effects of phase errors on 
the target localization MSE. First, it is apparent that the HCRB may be expressed as the sum of the CRB with no 
phase error and a term dependent on the statistics of the phase errors. This term is a function of the sensors location 
with respect to the target, through the matrices B m , and the system parameters (SNR, phase errors variance a\ and 
the number of mismatched transmitting and receiving radars) through the coefficients [i m . The manner in which the 
number of radars, their spread and the phase synchronization error variance affect the performance is not readily 
understood from (24). For this reason, numerical examples are employed in the next section to gain some insight 
into the relationships between system parameters and performance degradation. 



IV. Numerical Analysis 

We have evaluated the HCRB expression given in (20) numerically using the following example: M = 11, N = 9 
and a\ = [0,0.0001,0.0005,0.001,0.005,0.01,0.05], where a\ is expressed in (rad 2 ). The HCRB (a; , y Q ) is 
drawn in Figure 1. As a\ increases beyond a specific value, the additional CRB term ACRB dominates the 
performance and the curve. For high phase error levels, the performance degradation starts at lower SNRs. For 
small phase errors, localization accuracy is not undermined by the phase mismatch, and the HCRB (x ,y ) curve 
follows the CRB (x ,y ) closely. 
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case3-HCRB-10x.emf 



Fig. 1. HCRB for M=ll and N=9. The blue line represent the CRB value with no phase errors. c\ range from to 0.05. 



For a given system, the tolerated [u%\ max may be determined by solving ACRB ( [cr^] max ) ^ CRB q (x q , y a ). 
This value can serve as a design goal in the system phase calibration. For a given phase synchronization error 
variance u\, the expression ACRB (c^) gives the localization accuracy penalty. 

V. Conclusions 

MIMO radar with coherent processing exploits the signal phase measured at the receive antennas to generate 
high resolution target location estimation. To take advantage of this scheme, full phase synchronization is required 
among all participating radars. In practice, inevitable phase synchronization errors reflect on the system localization 
performance. In this paper, a closed-form expression of the HCRB of target localization has been derived, capturing 
the impact of the phase synchronization errors on the achievable target localization accuracy. In particularly it has 
been shown that the HCRB can be expressed as a sum of the CRB with no phase error and a term that represents 
the phase error penalty. The latter has been shown to be a function of the sensors geometry, SNR, and the number 
of transmitting and receiving radars in addition to the phase error MSE. As phase synchronization over distributed 
platform is a complex operation and phase errors are unavoidable, the HCRB offers valuable information at the 
system design level. For a given phase error MSE, the HCRB may be used to derive the attainable target localization 
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accuracy. Otherwise, for a given system performance goal on localization accuracy, the HCRB provides with an 
upper bound on the necessary phase error MSE values. 



Appendix I 
Derivation of the J(D) matrix 

In this appendix, we develop the elements of the matrix (k), i.e. [J^ (k)^ . = —E Kr \ Knr j 
based on the conditional pdf in (9). The diagonal submatrix R T is derived as follows: 

' d 2 In [p (r|«)) 



E r 



d 2 ln(p(r|K)) 



-E Kr \ Knr < E r \ K 



ol Re { W 



drikdr^k' 
s k (t-T lk )s* k , (t-T' tk ) 



and 



R t = 8tt 2 (/ c 2 + /3 2) snrlQxQ7 

where snr = |^| 2 /cr 2 , and the following notation is used: 

i = [(l-l)*M + k] and j = [(£' — l)*M + k'], 
k,k' = l,..,M, £, 
The elements of the matrix are given by 



1,..,N, 



_ 2MNsm _ 
t^VilA ~ 7772 ~ 1-^2,2 ' 



and 



[s tf ] 1)2 = = [S] 2)1 , 

and the elements of the matrix Sa are given by 

NImxm (H T ) 



Sa = 2 snr 
The off-diagonal submatrices are as follows: 



MxN 



(11^) 



NxM 



xN 



LxL 



F tAqxL = 47r/ c snr 

n(f)= [ o 



ImlQxl 


~^R 1qx 


Imxm 


n(i) " 


Imxm 


U(N) 



Nx(£-1) J-iVxl u Nx(N-£-l) 



QxL 
NxN 



and 



2 snr 



■pi 



d lm NlJ xM d lm MlJ xN 
-# Re NlJ xM -d Rc MlJ xN 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 
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